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Abstract 

Background: Multiple Sclerosis (MS) is considered a T-cell-mediated autoimmune disease with a prototypical 
oscillatory behavior, as evidenced by the presence of clinical relapses. Understanding the dynamics of immune cells 
governing the course of MS, therefore, has many implications for immunotherapy. Here, we used flow cytometry to 
analyze the time-dependent behavior of antigen-specific effector (T eff ) and regulatory (T reg ) T cells and microglia in 
mice model of MS, Experimental Autoimmune Encephalomyelitis (EAE), and compared the observations with a 
mathematical cross-regulation model of T-cell dynamics in autoimmune disease. 

Results: We found that T eff and T reg cells specific to myelin olygodendrocyte glycoprotein (MOG) developed 
coupled oscillatory dynamics with a 4- to 5-day period and decreasing amplitude that was always higher for the T eff 
populations, in agreement with the mathematical model. Microglia activation followed the oscillations of MOG- 
specific T eff cells in the secondary lymphoid organs, but they were activated before MOG-specific T-cell peaks in the 
GNS. Finally, we assessed the role of B-cell depletion induced by anti-CD20 therapy in the dynamics of T cells in an 
EAE model with more severe disease after therapy. We observed that B-cell depletion decreases T eff expansion, 
although its oscillatory behavior persists. However, the effect of B cell depletion was more significant in the T reg 
population within the CNS, which matched with activation of microglia and worsening of the disease. 
Mathematical modeling of T-cell cross-regulation after anti-GD20 therapy suggests that B-cell depletion may 
influence the dynamics of T cells by fine-tuning their activation. 

Conclusions: The oscillatory dynamics of T-cells have an intrinsic origin in the physiological regulation of the 
adaptive immune response, which influences both disease phenotype and response to immunotherapy. 

Keywords: T cells, Effector, Regulatory, B cells, Dynamics, Autoimmunity, Multiple sclerosis, Systems biology, 
Immunotherapy, Anti-CD20 



Background within the central nervous system (CNS), which pro- 
Relapsing-remitting oscillatory behavior is a hallmark duces demyelination and axonal damage and impairs 
of autoimmune diseases such as Multiple Sclerosis neural conduction, leading to clinical symptoms. Acute 
(MS) [1,2]. At the clinical level, the presence of re- inflammatory lesions are revealed by magnetic reson- 
lapses defines the subtypes of MS. The clinical relapse ance imaging (MRI) as contrast-enhancing lesions, 
rate is about 0.5-1 per year, with around a month of The number of clinically defined relapses represents 
duration and self-resolution. At the pathological level, only one in every 5-10 contrast-enhancing lesions, as 
a relapse is the result of an acute inflammatory process has been shown by MRI studies [3], because many 
brain regions are silent and brain plasticity copes 
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onset, although one third of the relapses are preceded 
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by common infections or stressful events [4,5]. Under- 
standing the biological basis of relapses in MS may 
have implications for immunotherapy. 

In autoimmune diseases, the balance between the pro- 
inflammatory response and regulatory factors seems to 
be disrupted. As such, antigen specific effector T cells 
(T e ff) are over-activated [6], whereas the function of regu- 
latory lymphocytes (T reg ) is altered [7-9], indicating im- 
pairment of peripheral tolerance. It could be expected that 
this impaired immune response would produce a chronic 
inflammatory process leading a progressive clinical course. 
However, in MS and other autoimmune diseases the pre- 
dominant behavior is fluctuating, with periods of T-cell ac- 
tivation and tissue damage followed by deactivation and 
repair. Although this behavior could be partially accounted 
for by external triggers (e.g. viral infections, stress, changes 
in light exposure and vitamin D levels), its robustness 
among pathologies, resilience in front of immunotherapies 
and the lack of closely related triggers suggest that this 
may not be the only explanation. We have postulated that 
the oscillatory behavior of the autoimmune response arises 
from an underlying periodic dynamics that is intrinsic to 
the design of the immune system, with several control 
mechanisms providing negative feedbacks that tightly con- 
trol T-cell activation [10]. Mathematical models predict 
that cross-regulation of T eff and T reg cells generates a 
stable oscillatory dynamic of both populations that main- 
tains homeostasis in health, but which also promotes 
relapsing-remitting flares under autoimmune conditions 
[10-12]. The generation of an autoimmune response 
seems to be based in the presence of susceptibility factors 
such as genetic polymorphisms, levels of modulators of 
the immune response such as vitamin D, or previous in- 
fections that would change the thresholds of T-cell activa- 
tion and differentiation, defining the autoimmune regime. 
However, the presence of environmental factors (e.g. vita- 
min D levels, acute infections or stress) or genetic poly- 
morphism have not been convincingly associated with this 
oscillatory behavior, supporting the concept that such os- 
cillations are inherent to the organization of the adaptive 
immune system. 

B-cell depletion has recently emerged as a powerful 
immunotherapy for treating autoimmune diseases such 
as Rheumatoid Arthritis (RA), Lupus or MS [13,14]. Al- 
though originally conceived for targeting the role of anti- 
bodies in the pathogenesis of autoimmune diseases, the 
quick and profound response observed in clinical trials 
in patients with MS and RA suggest that instead of de- 
creasing antibody production, B-cell depletion may work 
by modulating the pro-inflammatory environment cre- 
ated by B-cells, the antigen presentation function they 
provide to T-cells, or the function of regulatory B-cells 
[15,16]. Moreover, anti-CD20 therapy ameliorates the 
condition of the animal model of MS, Experimental 



Autoimmune Encephalitis (EAE), preventing the pro- 
inflammatory environment (Thl/Thl7 responses) pro- 
moted by B-cells [17]. However, in the EAE model in 
which B cells do not play an encepalithogenic role, B-cell 
depletion decreases T reg and regulatory B-cells, worsening 
the disease course [17]. Understanding the effects of B-cell 
depletion in T-cell activation would improve our under- 
standing of the mechanisms governing peripheral toler- 
ance and its role in autoimmune diseases. 

Here we validate at the experimental level the hypoth- 
esis that the relapsing-remitting behavior of autoimmune 
diseases is highly dependent on the cross-regulation of 
the lymphocyte populations, demonstrating that the 
underlying oscillatory behavior is intrinsic to the design 
of the immune system. We have analyzed the dynamics 
of Myelin Oligodrendrocyte (MOG)-specific T eff and 
T reg cells in secondary lymphoid organs (spleen) and in 
the target tissue (CNS), as well as the activation of 
microglia in the animal model of MS, the Experimental 
Autoimmune Encephalitis (EAE). Moreover, we used our 
experimental and mathematical model to obtain insights 
into the role of anti-CD20 therapy in the regulation of 
T-cell dynamics and the generation of autoimmune 
response. 

Results 

Antigen specific T cell subpopulations and microglia 
display coupled and interlinked oscillatory dynamics 
during the autoimmune attack in the brain 

In order to assess the presence of oscillation in T cell 
populations, their relationship and their behavior during 
the autoimmune response, we measured the frequency 
of both total and active MOG-specific T eff and T reg cell 
populations (compared with the non-specific T cells and 
the overall CD4+ cell population) using MOG-specific 
MHC class II tetramers during the course of EAE in the 
C57BL/6 mice immunized with MOG (Figure 1). T cells 
were measured in the spleen and CNS daily in series of 
animals (3 animals per day) immunized at the same time. 
As predicted by the cross-regulation model of T cell dy- 
namics [10-12], we observed that the activation of MOG- 
specific T eff and T reg cells in peripheral lymphoid organs 
(spleen) follows an oscillatory behavior (Figure 2A). Dur- 
ing the 30 days after immunization, we observed a consist- 
ent average time between peaks (period of oscillation) at 
around 4-5 days (Table 1), with coordinated dynamics be- 
tween the two T-cell populations. Also, the oscillation 
amplitude (maximum percentage of antigen- specific cells 
with respect to total CD4+ cells) of MOG-specific T eff 
cells was always significantly larger in average (between 2 
and 6 fold) than the one for MOG-specific T reg cells 
(Table 1). Additionally, the amplitude was always larger 
for the first two peaks, coinciding with the onset of clinical 
disease, and was attenuated along the course of the 
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Figure 1 (See legend on next page.) 
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(See figure on previous page.) 

Figure 1 MOG-specific T cells and microglia subpopulations analyzed by flow cytometry in mice with EAE. Leukocytes were obtained 
from the spleen and the CNS from MOG 35 _55 immunized C57BL/6 mice at different stages of EAE as indicated. A) T eff and T reg lymphocytes are 
shown. Numbers indicate the percentage of cells within the CD4+ cell gate. The histogram on the right shows the expression of CD69 on the 
surface of T eff and T reg cell subsets in the spleen. Numbers indicate the percentage of CD69+ cells. B) CD4+ cells stained with MOG 3 5_5 5 /IA b 
tetramer. MOG-specific T eff and T reg cells are shown. C) Dot-plots showing CD45 + CD1 1 b + MCH-ll(IA b ) + microglia population. D) Antigen-specific 
T eff and T reg oscillate over the course of EAE. Numbers indicate the percentage within the CD4+ cell gate. MOG-specific T cells from spleen (a, b) 
and CNS (c, d) in the onset, peak and recovery stages of EAE are shown. Animals were not treated (a, c) or treated with anti-CD20 antibody 
before immunization (b, d). 



disease. We observed that both T cell populations were 
phase-locked, with differences between the time peaks be- 
ing always smaller than 1 day, which indicates that the os- 
cillations of T e ff and T reg are tightly controlled by 
common mechanisms. In summary, after immunization, 
antigen-specific T e ff and T reg populations display an oscil- 
latory behavior. The robustness of the period observed 
also suggests the presence of a common regulatory mech- 
anism, supporting the cross-regulation model. 

We also measured the frequency of MOG-specific T e ff 
and T reg cells in the CNS of animals suffering EAE from 
the time at which cells migrate to the brain (from day 
9). Because activated T-cells migrate to the target tissue 
after T-cell activation, we predicted a similar oscillatory 
behavior in the CNS by infiltrating T cells unless the 
blood-brain barrier (BBB) was buffering these oscilla- 
tions. Accordingly, we observed the presence of MOG- 
specific T ef f and T reg cells in the CNS at various times, 
again showing a coupled oscillatory dynamics. Ampli- 
tudes of T-cell incursions were around 3% and 4% of the 
total population of non-specific T-cells for T e ff and T reg , 
respectively (Figure 2B, Table 1) and always significantly 
larger than the one observed in the spleen, indicating 
the accumulation of antigen-specific cells in the target 
tissue. This is in agreement with previous observations 
in the same animal model showing that T-cell invasion 
of the CNS starts 2-3 days after T-cell activation in the 
periphery (day 9-10 post-immunization) and extends 
along the course of the disease (Figure 2C) [18]. 

Considering that microglia are activated in response to 
CNS damage or by immunological mediators, we pre- 
dicted a similar oscillatory behavior of activated microglia 
(CD45 low CDllB + MHC-IA b ) population following CNS in- 
vasion by MOG -T-cells. Accordingly, we observed the 
presence of oscillations in microglia activation after EAE 
induction (Figure 2D). Activated microglia shows two 
peaks of activation with a period of 8 days and a max- 
imum amplitude of 60% (Table 1). Microglia remains acti- 
vated for 6 days in both activation events, indicating an 
intrinsic activation period for microglia as well, and differ- 
ent for the one observed for T cells (Figure 2D). However, 
microglia activation started at day 10, preceding the peak 
of MOG-specific T cells within the CNS, and exhibits a 
first maximum at day 14, one day before the following 



peak of MOG- T e ff into the CNS. At same time, the 
maximum of clinical relapse was registered on day 13 to 
14. The clinical score dropped from day 14 to day 21 
(Figure 2C), coinciding with the microglia deactivation 
to baseline levels. The second microglia activation peak 
spans from days 20 to 26, again preceding the following 
peak of activated MOG- T e ff cells in the CNS (Figure 2B). 
Overall, our data indicates activation of microglia preced- 
ing or concomitant with the invasion of CNS by MOG 
specific T-cells. In order to quantify the dependence of 
microglia activation on CNS infiltration by MOG-specific 
T-cells, we plotted the distribution of phases within a 
microglia cycle for which the average of MOG-specific 
T-cell incursion is observed in the CNS, observing again 
that peaks of microglia activation precede the infiltration 
of MOG-specific T cell to the CNS (Figure 3). Together, 
these findings indicate that microglia activation also os- 
cillates during brain autoimmunity, following MOG- 
specific T-cell activation in the periphery but preceding 
the incursion of those cells into the CNS. This result 
suggests the need for additional signals for microglia ac- 
tivation preceding antigen-specific T cell infiltration of 
the CNS and/or a role for microglia in the recruitment 
of MOG-specific T-cell within the CNS. 

Mathematical modeling identifies T reg activation as a 
major element governing autoimmune susceptibility 

In order to interpret the experimental data described 
above, we make use of a cross-regulation model of T-cell 
activation previously described (Figure 4A) [10] after 
adapting it to reproduce mice data (Additional file 1: 
Table SI). The model is based on the hypothesis of the 
strong influence of a negative feedback between T e ff (E) 
and T reg (R) cells in the adaptive immune response, using 
a predator-prey system (see equations in Figure 4A). Ac- 
tive T e ff would produce signals to evoke expansion of T reg 
population; meanwhile, the increase of T reg activity would 
mean a suppressive effect on the levels of active T e ff. The 
model also includes thymic resting/naive T cells (R r and 
E p Figure 4A) arising from a discrete stochastic input of 
cells in the immune system, mimicking the apparently 
random nature of the pulse trains of resting cells (which 
can become activated by antigen presentation). We identi- 
fied two key characteristics of the cross-regulation model 
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Figure 2 Dynamics of activated MOG-T eff , MOG-Treg cells and microglia in secondary lymphoid organs and CNS during the course of EAE. 

A) Percentage of MOG-specific T eff (red) or T reg (blue) in the spleen after immunization. The upper graph shows the mean of the measurements from 
3 animals each day for both populations (T eff and T reg ) together. The bottom panels show the distribution of all measurements (raw data) from 
3 animals per day for each subpopulation as dots, with the mean in black; B) percentage of MOG-specific T eff (red) orT reg (blue) in the CNS from day 9 
after immunization. The upper graph shows the mean of the measurements from 3 animals each day, while the bottom panels show the distribution 
of all measurements for each subpopulation as dots, with the mean in black; C) Clinical score of animals suffering EAE after immunization (mean + SD); 
D) Percentage of activated microglia in the CNS from day 9 after immunization. 
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Table 1 Statistics of period and amplitude of the oscillation of MOG -specific T cells and microglia in EAE 



Period Amplitude Period Amplitude 

(days) (%) (days) (%) 
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0.06 ± 0.02 * 


CNS 


MOG-T eff 


4± 1 
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CNS 


MOG-T eff 


3.0 ± 0.4 


2.4 ± 1 .4 




MOG-T reg 


2.3 ± 0.9 


4.0 ± 0.9 *** 






MOG-T reg 


3.3 ±0.7 


0.5 ± 0.2 *** 




Microglia 
(activated) 


8 

Period 
(days) 


60 ±7 

Amplitude 
(cell num) 






Microglia 
(activated) 


Period 
(days) 


57±9 

Amplitude 
(cell num) 


Model EAE Spleen 


MOG-T eff 


5.2 ±0.2 * 


2254 ±128 *** 


Model EAE anti-CD20 


Spleen 


MOG-T eff 


6.5 ± 0.5 * 


816 ± 18 *** 




MOG-T reg 


4.9 ± 0.2 


552 ± 26 *** 


therapy 




MOG-T reg 


5.5 ±0.4 


429 ± 1 0 *** 


Period and amplitude of the peaks of active MOG-T e ff , 


MOG-T r eg and activated microglia in secondary lymphoid organs (spleen) and the CNS in animals 



developing EAE (as described in Figure 2), in model simulations (as described in Figure 4) and animals with EAE and treated with anti-CD20 therapy (as described 
in Figure 6). Amplitudes of MOG-specific T cells and activated microglia are shown as percentages to the CD4 + population and total microglia population, 
respectively; amplitudes in model simulations are expressed as cell number values. Data is shown as the mean ± SEM, statistical significance are between EAE 
versus EAE anti-CD20 therapy; sig. (*) p < 0.05, (***) p < 0.0001 (paired T test). 



that clearly account for the difference between the healthy 
and EAE behavior, namely the average amplitude of the 
peaks and the average periodicity. Using those features, we 
established that the rate of T reg activation (a reg ) and the 
number of T eff cells required for half-maximal activation 
of T reg cells (Keff ) were found to be critical for reproducing 
the oscillatory autoimmune behavior data versus the 
healthy configuration parameters. In particular, using the 
same value for the rate of antigen presentation in both 
sub-populations (5), we observed oscillatory behavior after 
a decrease of a reg (down to 1.05 /day; range: 0.75 to 1.05) 
and an increase of K e ff (up to 2,500 cells) (Figure 4 B, C). 
Thus our data and model sensitivity support the concept 
that the observed oscillatory autoimmune response could 
be explained by a decrease in the strength of the T efr T reg 
cross-regulation, a combination of increasing the T e ff 
threshold for activating T reg cells and reducing the num- 
ber of activated T reg cells. 

In agreement with our experimental data, the mathem- 
atical simulations showed that the coordinated oscillations 



of the active populations of T eff and T reg cells have ampli- 
tudes for T eff at least 4 times larger in average than T reg 
(Table 1). Also, in order to compare the adapted model to 
mice conditions with our experimental data, we tested the 
average distribution of standard deviations against the 
means of the oscillation amplitudes for three simulations 
(1 point per day) for 30 days. Means of relative deviations 
measured in the experimental data was 0.62 ± 0.01 for T eff 
cells and 0.40 ± 0.02 for T reg cells, while in model simula- 
tions it was 0.61 ± 0.01 for T eff and 0.30 ± 0.01 for T reg 
cells. These results fit qualitatively with the experimental 
observations in the EAE dataset. 

The model proposes a predator-prey system in which 
T reg and T eff play the active-subsets of predator and prey 
populations, respectively. We performed a phase-space 
representation of the mathematical model trajectories in 
the plane for T eff -T reg after fitting to experimental EAE 
data (Figure 4D), showing that the healthy-state dynamics 
are controlled by T reg (x-axis), while the autoimmune- 
state dynamics are mainly governed by T eff (y-axis). In 
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Figure 3 Analysis of the time dependence between microglia activation and MOG-specific T-cell infiltration of the CNS. Distribution of 
the average of incursion of MOG-specific T cells into the CNS (Jeff and Treg) versus microglia phase cycle. 
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Model: ODE system for T cells populations 
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summary, the combination of experimental data with 
theoretical modeling supports the importance of cross- 
regulation of T-cell populations in the control of the im- 
mune response, and points to T reg cells as a critical popu- 
lation in order to understand dynamics of autoimmune 
diseases. 

B cell depletion modulates T-cell subpopulation dynamics 
differentially, influencing the outcome of EAE 

Recently, it has been shown that B-cell depletion by anti- 
CD20 antibodies has a strong effect in suppressing the 
relapsing-remitting course of autoimmune diseases such 
as RA and MS [13,14], modulating the frequency of T-cell 
subpopulations [15]. In order to analyze the influence of 
B-cells on T-cell dynamics, we analyzed the frequency of 
antigen-specific T eff and T reg subpopulations in mice dur- 
ing the course of EAE after B-cell depletion. We made use 



of C57BL/6 mice immunized with MOG peptide, given 
that in this model B-cells are not encephalitogenic, but ra- 
ther have been shown to play a role in T cell activation 
[17]. As expected, we observed that anti-CD20 therapy in- 
duced a rapid, almost complete depletion of B-cells in the 
peripheral lymphoid organs (Additional file 1: Figure SI). 
Regarding the dynamics of T-cells, we observed that after 
B-cell depletion, MOG-specific T eff and T reg subpopula- 
tions maintained their oscillatory presence in secondary 
lymphoid organs, although in this case the large activation 
peaks of T eff that appeared in non-treated animals suffer- 
ing EAE were significantly decreased in amplitude and 
appeared at earlier times (Figure 5A). However, similarly 
to untreated animals suffering EAE, the peaks of T eff and 
T reg were phase-locked. Interestingly, in the CNS T eff and 
T reg incursions were present as well, but with a significant 
decrease in the amplitude of T reg cells compared to T eff 
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cells (Figure 5B, Table 1). Surprisingly, even if the fre- 
quency of both MOG T-cell subpopulations were de- 
creased, microglia were activated to a similar extent as in 
untreated animals with EAE (-60% in the first peak with a 
maximum at day 16), followed by oscillations without a 
clear period or phase, and did not fall to baseline levels 
(Figure 5D), suggesting a transition to chronic activation 
and contributing to the worsening of the disease. At the 
same time, the maximum clinical score was delayed com- 
pared with untreated animals with EAE (day 17 in Fig- 
ure 5C compared to day 12 in Figure 2C) and disease was 
more severe, as previously described in this model [17]. 
Both groups showed statistically significant differences by 
paired-t test (with p < 0.05), specifically at days 11 (early 
phase of the disease episode), 17 (maximum score of 
treated EAE group, Figure 5C), at 21 (late phase of the dis- 
ease episode). 

The presence of an imbalance in the CNS of the peaks 
of antigen-specific T eff compared to T reg and the persist- 
ence of microglia activation may account for the impair- 
ment of EAE in this model despite anti-CD20 therapy. 
Therefore, our results suggest that B-cell depletion mod- 
ulates T-cell activation in peripheral immune organs 
with more profound effects in T reg activation (Table 1, *, 
p < 0.05). Downregulation of T reg was mainly revealed in 
the CNS tissue (Table 1, p < 0.0001), and was associated 
with persistent microglia activation and (as mentioned 
before) with the higher clinical severity in the pre-treated 
EAE group. However, the fact that the oscillatory dynamics 
of MOG-specific T-cells remain suggests that B-cell deple- 
tion works by tuning T-cell activation thresholds, rather 
than reverting the autoimmune process itself. 

Sensitivity analysis of the cross-regulation model support 
the critical role of B-cell depletion on T reg activity 

In order to analyze the influence of B-cell depletion on 
the dynamics of T-cell subpopulations, we performed a 
sensitivity analysis of our mathematical cross -regulation 
model with the aim of identifying the model parameters 
governing the dynamics induced by this therapy in T 
cells. Although our model does not include B- cells, by 
fitting the model with the experimental data we were 
able to analyze the resulting dynamics of both T-cell 
populations after anti-CD20 therapy and obtain insights 
into a possible mechanism for explaining how B-cells in- 
fluence T-cell dynamics. We found that the dynamics of 
the antigen-specific T-cell subpopulation after anti-CD20 
therapy was reproduced by reducing the K eff threshold 
below the healthy standard (< 1,000 cells; e.g. 850 cells), in- 
dependently of the a reg parameter (a reg values were in the 
low range of autoimmune regime; Figure 6). In other 
words, our mathematical model suggests that B-cell de- 
pletion therapy may influence the autoimmune process 
by preventing uncontrolled activation of T eff without 



strengthening T reg activation. Experimental trajectories 
through the phase-space were dominated by the T eff cells 
(y-axis, Figure 6C) in EAE (in blue) and in anti-CD20 
therapy (in red). However, in the anti-CD20 treated 
group, we did not observe the high excursions character- 
istic of EAE. This behavior can be easily modeled by 
performing simulations for 30 days, which highlights the 
fact that B cell depletion, although suppressed the T eff ex- 
cursions, did not recover the T reg excursions observed in 
the healthy cohort (Figure 6D). 

We performed simulations of the changes in the survival/ 
death rate of T cell populations in order to predict its 
effects in the oscillatory dynamics of T cells. We observed 
that a 2-fold decrease or increase in the cell death rate of 
T eff or T reg decreases the period by at least one day for 
active T e ff cells and two days for T reg . Although peak 
period (period of time between two consecutive max- 
imums of activation) is not a well-known characteristic of 
T cell activation, it would be a good marker to confirm the 
oscillatory nature of this process. In summary, our simula- 
tions suggest that B-cell depletion modulates autoimmune 
diseases mainly by affecting the interaction between T reg 
and T eff , and in particular strengthening the T eff -T reg 
cross-regulation by lowering the T reg activation threshold. 

Discussion 

In this study we have assessed the dynamics of antigen- 
specific T-cell subpopulations after immunization and 
related them to the outcome of brain autoimmune disease. 
We observed that, after immunization, antigen-specific ef- 
fector and regulatory T-cells are activated and oscillate in 
a coupled manner with an intrinsic period. The combin- 
ation of experimental and computational analysis supports 
the concept that relapses in autoimmune diseases are an 
intrinsic property of the regulation of the immune system, 
that can be modulated by environmental factors [10]. The 
fact that oscillations of both populations were matched 
day-to-day (phase-locked), indicates the importance of the 
synchronous activation of effector and regulatory adaptive 
immune response for the outcome of the immune 
response [19]. Previous models of T-cell activation suggest 
that regulatory T-cell activation follows effector T-cell ac- 
tivation with some delay [20-22]. However, the fact that 
our experimental data shows that both regulatory and ef- 
fector T-cells become activated the same day after 
immunization points the importance of triggering the 
regulatory T-cell response at the same time and in parallel 
than the effector response. After activation, the negative 
feedback between both populations could maintain their 
oscillation for long period of time, although removal of the 
antigens and other suppressive factors cancels the oscilla- 
tions in the long term. The analysis of the T cell cross- 
regulation model after fitting with the experimental data 
supports a key role of regulatory T-cell dynamics in 
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defining the outcome of the autoimmune response [22,23]. 
Our model could be further validated by future experimen- 
tal studies testing, for example, different scenarios of T-cell 
viability assays (involving e.g. cellular pathways implicated 
in the control of cell death and survival, such as Jak/Stat, 
Ras/MAPK, PI3K/AKT, and interleukin signaling). 

During the course of EAE, activated CD4+ and CD8+ 
cells enter the brain, which is associated with microglia 
activation that mediates demyelination and axonal loss 
[18]. In the EAE model, regulatory T-cells have a critical 
role in controlling the activation of encephalitogenic 
T-cells, but their mechanism of action is complex. 
Regulatory T cells are dependent on the local (CNS) con- 
ditions that regulate their activation and suppressive activ- 
ity [24,25]. However, the functional outcome seems to be 
dependent on balancing the suppression of the activity 
and proliferation of effector T-cells with the promotion of 
the survival and activity of encephalitogenic T-cells, 



causing their accumulation and distorting the regulatory 
and effector T-cell balance [26]. 

Another interesting observation from our results is 
that the activation of the innate immune system, namely 
microglia activation, also follows the oscillatory dynam- 
ics of the adaptive immune response with a well-defined 
phase. Interestingly, we observed that microglia become 
activated closely after activation of T-cells in the periph- 
ery, but before the infiltration of antigen-specific T-cells 
into the CNS. This results points to a model in which 
the autoimmune response organized in secondary lym- 
phoid organs or in the circumventricular organs of the 
brain delivers signals (most likely soluble factors) to the 
brain parenchyma, activating microglia. Time series ana- 
lysis of gene expression changes along the course of 
EAE have shown activation of genes related with the 
innate immune response before the presence of inflam- 
matory infiltrates, supporting the role of microglia 
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activation in the early stages of the disease, even preced- 
ing an antigen-specific T-cell invasion of the CNS [27]. 
Recent pathological studies in the brain of patients with 
MS as well in EAE models suggest that activated effector 
T-cells accumulate in subaracnoid space where they can 
be restimulated, releasing pro-inflammatory signals that 
activates microglia and contributes to the opening of the 
blood-brain barrier and subsequent infiltration into the 
brain by T-cells of the CNS parenchyma [28-30]. 

Amplitudes of MOG T-cell activation peaks were clearly 
larger in the early stages of the disease following immu- 
nization and then decreased. This could be explained by 
several factors. After immunization, the availability of anti- 
gens at site of immunization to be presented to T-cells by 
antigen presenting cells is significantly higher compared 
with the antigen concentration within the CNS in 
consecutive restimulations. Also, immunization exerts a 
synchronization effect that makes it easier to detect the 
activation of T-cells right after the animal is immunized, 
in comparison with later times at which individual clones 
become increasingly desynchronized, complicating the de- 
tection of T-cell expansion. Another factor may be the ef- 
fect of epitope spreading, which involves the recruitment 
of other antigen specificities that were not measured here. 
Also, anti-inflammatory signals from the tissue, such as 
IL-10 would contribute to shutdown the immune response 
and therefore dampening of the peaks of activated T-cells. 

B-cells seem to have a complex role in brain auto- 
immunity. In addition to their role as antigen producing 
cells, B-cells have been shown to contribute to the 
pathogenesis of autoimmune disease via the production 
of a pro-inflammatory factors such as IL-6, the ac- 
tivation of pro-inflammatory Thl and Thl7 cells and 
monocytes, and the inhibition of T-cell activation by 
regulatory B (B10) cells [15,17,31-34]. The role of B-cell 
depletion in the activation of antigen specific T-cells in 
EAE is complex as B-cells do not directly impair T-cell 
function, but rather alter the balance between effector 
and regulatory function of these cells [35-37]. Our re- 
sults indicate that in a model in which anti-CD20 ther- 
apy induces a more severe disease course, even if B-cell 
depletion decreases the activation of both effector and 
regulatory T cells, the disease outcome seems to be 
dependent, at least in part, on the balance of these two 
populations. We show that the activation of MOG- 
specific T-cells happens earlier and faster in anti-CD20 
treated animals than in non- B-cell depleted animals. 
This is agreement with the view that effector B-cells 
modulate antigen presentation to T cells and regulatory 
B-cells modulate T cell expansion [37,38]. In the CNS of 
B-cell depleted animals, we observed a more significant 
effect on regulatory than effector response and a subse- 
quent enhancement of microglia activation, which corre- 
lates with more severe CNS tissue damage. 



Our results have several implications for immunother- 
apy. First, immunomodulatory therapies, such as anti- 
CD20 therapy, do not revert the pathogenic response in 
autoimmune diseases, but rather maintain them in a dy- 
namic state that is less deleterious. For this reason, ces- 
sation of the therapy would allow to the immune system 
to come back to its original autoimmune dynamics and 
in some cases, even worsen the disease due to a rebound 
effect. This has been well documented in cessation of 
natalizumab therapy in humans, in which a significant 
proportion of patients develop a rebound relapse a few 
months after the last natalizumab infusion. Second, al- 
though regulatory T-cells are an ideal target for immuno- 
therapy - e.g. increasing their function for treating 
autoimmune disease or decreasing their function for can- 
cer immunotherapy - the outcome of the disease will be 
dependent on the dynamics of all populations implicated 
and the timing in which the therapy is started. For this 
reason, the ideal therapy at this level would involve fine- 
tuning the long-term dynamics of antigen-specific T-cells 
to maintain them at a level closer to the healthy state. 

Our study has several limitations. First, we have quan- 
tified antigen-specific T-cells using MHC class II tetra- 
mers by flow cytometry, which has a resolution for 
identifying 1 in every 10 4 target cells in the tissue from a 
standard acquisition of 2-5 x 10 5 leukocytes [39]. More- 
over, because we attempted to quantify MOG specific 
T cells within the context of the natural repertoire of a 
rodent, instead of using transgenic animals with single 
TCR specific for MOG (which may have simplified the 
identification of MOG specific T-cells), quantification was 
done in different animals and individual dynamics were 
extrapolated from the population analysis. For this rea- 
son, we analyzed results as proportions to the reference 
CD4+ cell population instead using absolute numbers. 
While the presence of noise created by technical limita- 
tions might have prevented from observing more cycles 
due to desynchronization of the T-cell activation be- 
tween animals, the oscillations reported here were robust 
enough to be captured with this technique. Also, we fo- 
cused in the immune response specific to the single anti- 
gen used for immunization (MOG), but did not assess 
the role of epitope spreading in the dynamics of auto- 
immune T cells nor in the dynamics of other lymphocyte 
subpopulations (e.g. CD8+ cells) with other specificities, 
which may have a role in the outcome of the auto- 
immune response. Finally, anti-CD20 produces profound 
but not complete B-cell depletion, and does not target 
only MOG-specific B-cells. For this reason, residual 
B-cells may have played a role in the worsening of 
EAE. Alternatively, anti-CD20 therapy may preferentially 
target some B-cell subpopulations such as IL-10- or IL-6- 
producing B cells, or cells at different stage of differenti- 
ation that may have a specific effect on T-cell dynamics. 



Martinez-Pasamar et al. BMC Systems Biology 2013, 7:34 
http://www.biomedcentral.eom/1752-0509/7/34 



Page 12 of 15 



However, the results from our study support previous 
models of the regulation of T- cell populations and their 
role in the immune response. 

Conclusion 

Our study supports the importance of the cross - 
regulation between lymphocyte populations and their 
contribution to the outcome of the immune response. 
Moreover, highlights the role of regulatory T-cells in 
driving the dynamics of the immune response. Even if 
these cells are a smaller population compared with naive 
and effector cells, their critical role in the regulation of 
immune response makes them central in the pathogen- 
esis of autoimmune diseases, and thus attractive targets 
for immunotherapy against these diseases, even in an in- 
direct way such as is the case in anti-CD20 therapy. 

Methods 

Animals and induction of EAE 

EAE was induced in C57BL/6 (I-A b ) mice by injecting 
them with an emulsion containing 300 \i% of MOG 35 _55 
peptide and 500 \ig of M. Tuberculosis extract in incom- 
plete Freund adjuvant subcutaneously into the flanks as 
described before [40]. Mice receive 0.2 ml of the emul- 
sion in the flank. In addition, the mice receive 500 ng of 
B. Pertussis toxin via intraperitoneal injection (i.p) in 
200 \i\ PBS on days 0 and 2. Clinical signs of EAE were 
assessed according to the following score: 0, no signs of 
disease; 0.5, partial loss of the tone in the tail; 1, loss of 
tone in the tail; 2, hind limb paresis; 3, hind limb paraly- 
sis; 4, tetraparesia; 5, tetraplegia; 6, moribund [6]. Mori- 
bund mice were given disease severity scores of 6 and 
euthanized. For each experiment, we made use of 3 ani- 
mals per day (or every other day for repetitions) for 30 
days, and the experiments were repeated twice. The 
study was approved by the ethical committee on animal 
research of the University of Barcelona. 

Tissue preparation and T-cell isolation 

Splenocytes were obtained from the spleen by digesting 
it with collagenase D (Roche) and Dnase I (Roche) at 37°C 
for 45 min. Mononuclear cells were isolated by passing 
the tissue through a cell strainer (70 |im) followed by a 
Ficoll (Sigma) gradient centrifugation. T cells from the 
CNS were obtained by collecting the forebrain, cerebellum 
and spinal cord. CNS tissue was cut into small pieces and 
digested with collagenase D (Roche) and Dnase I (Roche) 
at 37° C for 45 min. Mononuclear cells were isolated by 
passing the tissue through a cell strainer (70 \im) to obtain 
single cell suspensions. Leukocytes were isolated from the 
CNS by gradient centrifugation. Briefly, a Percoll (Sigma) 
gradient (70/37%) centrifugation was made and inter- 
phase between 70% and 37% phase was taken. Myelin in 
the upper layer was removed. Cells harvested from the 



gradient inter-phase and the upper-phase was washed in 
PBS and resuspended. 

Tetramers purification and cell staining 

MOG 3 5_5 5 /IA b tetramer construct was generously pro- 
vided by Prof. Vijay Kuchroo, from Harvard University, 
and purified as previously described [25]. Tetramers 
were incubated with PBS, 0.2% BSA, 0,1% sodium azide 
for three hours at 37°C at darkness. After washing, cells 
were stained with 7-AAD, (BD Pharmingen) and anti- 
bodies against CD4 (BD Pharmingen), CD62L (BD 
Pharmingen), CD25 (BD Pharmingen), CD69 (BD 
Pharmingen), and CD45 (BD Pharmingen). For microglia 
activation, cell were stained with anti-MHC class II (IA b ) 
(Abeam), CD lib (BD Pharmingen) and CD45 (BD 
Pharmingen). B-cell staining was performed using anti 
CD45R/B220 (BD Pharmingen) and anti-CD21 (BD 
Pharmingen) antibodies. Stained cells were analyzed on a 
FACSCanto machine (BD biosciences) and data analysis 
was performed with FACS Diva software. 

Lymphocyte and microglia subpopulations analysis 

Antigen specific T cells were characterized by being tetra- 
mer positive (IA b - MO G+)- MOG-specific T eff cells were 
gated as the CD45 + CD4 + CD25"CD69 + IA b - MO G + popula- 
tion [25,41-43] (Figure 1A). MOG-specific T reg cells were 
gated as the CD45 + CD4 + CD25 hi IA b - MOG + population 
[8,44,45] (Figure IB). We did not check Foxp3 expression 
on the T reg population because it requires cellular perme- 
abilization, which was not compatible with the tetramer 
staining. Nevertheless, the subset analyzed corresponds to 
T reg population as described before [25] . Also, we analyzed 
the expression of CD69 and CD62 since there is a sub- 
population of T reg cells that expressed CD69. Non- 
encephalitogenic T eff lymphocytes were characterized by 
being CD45 + CD4 + CD25TA b - MO G~ (Figure 1C). Non- 
encephalitogenic T reg lymphocytes were characterized by 
being CD45 + CD4 + CD25 hi IA b - MO G~ (Figure ID). CD69 
and CD62 expression was also analyzed in non- 
encephalitogenic T-cells. B-cells were analyzed after 
staining with B220 and CD21 antibodies [31,46]. Deple- 
tion of mature B-cells (B220 + CD21 + ) was examined in 
the spleen and CNS every day over the disease course in 
the group of mice that received anti-CD20 treatment 
[17]. Activated microglia was defined as the CD45 low 
CDllB + MHC-II-IA b+ population (Figure 1C) [47,48]. 
MOG-specific T cells and non-encephalitogenic T cells 
(either T eff or T reg ) are shown as percentages with re- 
spect to the CD4 + population (MOG-Teff /CD4 + ; MOG- 
Treg/CD4 + ; Teff/CD4 + ; Treg /CD4 + ) in order to adjust 
for differences in total cell count per day. B lymphocyte 
results are shown as percentages of total lymphocytes 
(B Lymphocytes/ Total Lymphocytes). Microglia activa- 
tion is shown as the IA b percentage of total microglia. 
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Anti-CD20 therapy 

Mouse anti-CD20 monoclonal antibody was generously 
provided by Andrew Chan from Genentech. To deplete 
B-lymphocytes in vivo, animals received anti-CD20 anti- 
body 200 ug in 250 ul PBS one week before immuniza- 
tion as described before [17,31]. 

Mathematical model, simulations and sensitivity analysis 

We make use of a model of active T eff -T reg cross - 
regulation developed previously to describe the dynam- 
ics of T-cells in humans [10], updating here the model 
parameters to reproduce experimental data from EAE 
studies in mice. The model is based in 4 differential 
equations describing the dynamics of antigen specific 
resting T eff (1), resting T reg (2), activated T eff (3), and ac- 
tivated T reg (4). 

dE 

-^ = I E -E r S-E r p + Erj (1) 



Statistical analysis 

All experiments were performed at least two times. The 
values are expressed as the means ± SEM. We use the 
Mann-Withney test to determine statistical significance 
between groups and Wilcoxon signed-rank test for dif- 
ferences in time series. Statistical analyses were per- 
formed using PASW 18.0 software (IBM) and GraphPad 
Prism 4 software. 

Additional file 



Additional file 1: Table SI. Parameters of the model and initial 
conditions. Figure SI. B cell counts after anti-CD20 therapy. Percentage 
of B cells in the spleen (L-B spleen) and the CNS (L-B CNS) from two 
immunized animal per day are plotted along the duration of the 
experiment (30 days). Figure S2. Predictions of the effects of changes in 
the survival/death rate of T cell populations in the oscillatory dynamics of 
T cells. Graphs shows simulations of the effect of a 2-fold decrease in the 
cell death rate of T eff and T reg in the period of T cell population 
oscillations. 
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dR r /_x MS: Multiple Sclerosis; EAE: Experimental Autoimmune Encephalomyelitis; 

~~jf = iR-KrO-KrP + Krj (2) j^. E ff ec tor T cells; T reg : Regulatory T cells; MOG: Myelin olygodendrocyte 

glycoprotein; MRI: Magnetic resonance imaging; RA: Rheumatoid Arthritis. 
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Specifically, parameter changes leading to time rescal- 
ing were performed in order to adapt the characteristic 
period of T-cell oscillations from the values observed ex- 
perimentally, with respect to previous theoretical values 
for humans [10] (Additional file 1: Table SI). In this new 
version of the model we also represent the pulse train 
input of naive T-cells from the thymus by a noise-driven 
excitable dynamics typical of some gene circuit architec- 
tures, as a mathematical tool to mimic continuous ran- 
dom activity pulses [49]. We simulated all stochastic 
differential equations of the model by means of a sto- 
chastic Runge-Kutta method of second order (Heun 
method) using Matlab software. 

We performed a sensitivity analysis of the model by 
increasing and decreasing the baseline parameters 
(Additional file 1: Table SI) by 20% and 50%, and 
collecting the statistics of the average amplitude and 
period of activated T e ff and T reg populations. This sensitiv- 
ity analysis was performed by running simulations over a 
time of 300 days in model units using the Matlab software. 
Simulations are initiated with no resting T e ff and T reg cells, 
and 1,000 and 200 activated T eff and T reg cells, respect- 
ively, as described before [10]. 
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